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' 1. Introduction 

This paper analyzes the finite horizon optimal stopping problem for an n-dimensional jump diffusion process X 
which is governed by the following stochastic differential equation: 

(1.1) dX t = b(X t -,t) dt + a(X t -,t) dW t + dj u 

in which W = {Wt]t > 0} is the d-dimensional standard Brownian motion and J = {Jut > 0} is a pure jump 
Levy process independent of the Brownian motion. This jump process J can be of finite/infinite activity with 
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finite/infinite variation. We denote the Levy measure of J as v (please refer to Section 2 for the definition of J 
and its properties). 

We investigate the problem of maximizing the discounted terminal reward g by optimally stopping the process 
X before a fixed time horizon T. The value function of this problem is defined as 

(1.2) u{x,t)= sup E*'* [e- rr g(X T j\ , 

T£7o,T-t 

in which 7o,t is the set of all stopping times valued between and t. A specific example of such an optimal stopping 
problem is the American option pricing problem, where X models the logarithm of the stock price processes and 
g represents the pay-off function. In [1], Ait-Sahalia and Jacod consider the model in (1.1) and find evidence of 
infinitely active jumps in stock prices. 

The function u satisfies, a variational inequality with a nonlocal integral term (see e.g. Chapter 3 of [6]). In 
the literature different solution concepts were studied. Pham showed in [27] that the value function of the optimal 
stopping problem for a controlled jump process is a viscosity solution of a variational inequality using the dynamic 
programming principle. Also see [2], [3] for more recent results in this direction. In [23], Lamberton and Mikou 
proved that the value function associated to the optimal stopping problem for Levy processes can be understood 
as the unique solution of the same variational inequality in the distributional sense. In general, without assuming 
that the diffusion coefficient is non-degenerate, the value function is not expected to be a smooth solution of this 
variational inequality. 

Using the non-degeneracy assumption, Sections 1-3 in Chapter 3 of [6] and in [16] developed the regularity 
results for the Cauchy problem and boundary value problems for second order partial integro-differential equations. 
They proved existence and uniqueness of solutions in both Sobolev and Holder spaces. Also see [24]. On the 
other hand, there are only limited results for obstacle problems cither finite or infinite activity jumps. Bcnsoussan 
and Lions showed in Theorem 4.4 pp.250 of [6] that the solution of a variational inequality on a bounded domain 
can be characterized as an element in a certain Sobolev space. Their regularity results are not enough to ensure 
the smooth- fit property. Later, these results were extended to problems on unbounded domains by [19] and [32], 
but these papers assumes that the state variable X follows a diffusions or a jump diffusions with finite activity 
jumps. [26], [31], and [4] further analyzed the variational inequality for jump diffusions with finite activity jumps 
using different techniques. More recently, [11], generalizing the results in [18] for the diffusion case, analyzed the 
regularity of the value function of infinite horizon impulse control problems for infinite activity but finite variation 
jumps. A regularity result which treats the obstacle problem on unbounded domains with infinite variation jumps 
has been missing in the literature. 

In this paper we extend the results of [6] and analyze the optimal stopping problem for Levy processes with 
infinite activity jumps on an unbounded domain. In our main result (Theorem 4.1), we show that the value function, 
which is the unique solution of a variational inequality, is an element of Wp'i oc (see Section 4.1 for the definition 
of this Sobolev space). This regularity result directly implies that the smooth fit property holds and the value 
function is smooth inside the continuation region. 

We use the penalty method to analyze the variational inequality. The infinite activity unbounded jumps bring 
technical difficulties in applying this method. In particular, the usual L p ~ estimates are no longer directly applicable. 
Our main technical contribution is the norm estimate in Theorem 5.1, which we prove using an interior estimation 
technique. 

In Section 3, we treat the finite variation jump case separately because the reduced form of the integral operator 
and the Lipschitz continuity of the value function enable us to use classical arguments to prove that the value 
function is C 2 ' 1 the continuation region. We present this section because contrasting it with the later sections 
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clearly demonstrates the technical difficulties of handling the infinite variation case due to the complicated form of 
the integral operator in this case. 

The rest of the paper is organized as follows. In Section 2, we introduce the variational inequality and recall 
two notions of generalized solutions studied in [27] and [23]. In Section 3, we discuss the finite variation jump case 
and analyze the regularity of value function in the continuation region. Section 4 is devoted to study the regularity 
when jumps may have infinite variation. The main result of this paper (Theorem 4.1) is proved in Section 5. A 
key estimate (Theorem 5.1) is proved in Section 6. Proofs of several auxiliary lemmas are given in the Appendix. 

2. Preliminaries 

2.1. A priori regularity of the value function. Let us first analyze the pure jump component J in (1.1). 
According to the Levy-Ito decomposition (see e.g. Theorem 19.2 in [28]), J can be decomposed as 

(2-1) J t = jf + )imJ?, 

in which 

(2.2) J?= f f yn(ds,dy), Jt=\ f yji(ds,dy), 

JO J\y\>l Jo Je<\y\<l 

represent large and small jumps respectively. Here is a Poisson random measure on R + x (R" \ {0}). Its mean 
measure is the Levy measure is, which is a positive Radon measure on K" \ {0} with a possible singularity at 0. 
Even with this possible singularity at 0, the measure v still satisfies 

(2.3) J (\y\ 2 M)v(dy) < +co. 

1/2 

Here, the norm | • | is the standard Euclidean norm: \y\ = (Y^^iiy 1 ) 2 ) ■ ^ n (2-2), fl{ds, dy) = fi(ds, dy) — ds v(dy) 
is the compensated Poisson measure. 

We assume that the drift and the volatility in (1.1) are bounded and Lipschitz continuous, i.e., there exists a 
positive constant Lb, a such that 

\b(x,t) - b(y,t)\ + \a(x,t) - <r(y,t)\ < L h . a \x - y\, Vx, y G R n , 

(HI) 

moreover, \b(x,t)\ and |cr(x, f)| are bounded on R™ x [0, T]. 

We name the solution of (1.1), with the initial condition Xq = x, as X x . Thanks to (HI), X x has the following 
norm estimates. 

Lemma 2.1. Let us assume b and a satisfy (HI). Then there exists a positive constant C such that 

(2.4) ~E\X*-X*\<C\x-y\, for any r G %,t with t < T and x,y eR™ . 
Moreover, if the Levy measure satisfies 

(H2) f \y\u(dy) <+oo, 

J \v\>i 

then we have 



(2.5) E\X?\<C, E\X* -x\ < Ct 1/2 , and E 



sup \Xf — x\ 

0<s<t 



<Ct 1 ' 2 . 



Remark 2.1. Similar estimates were given in Lemma 3.1 of [27] under a slightly stronger assumption on the large 
jumps: Jiy^i \y\ 2 v {dy) < +oo. Using the equivalence between the norm \y\ and the norm Y^i=i \y l \> one cou ld 
prove Lemma 2.1 under assumption (H2). We give its proof in Appendix A. 
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Let us assume that the terminal reward g : R™ — > R is a bounded and Lipschitz continuous function, i.e., there 
exist positive constants K and L such that 

(H3) < g{x) < K and 

(H4) \g(x)-g(y)\<L\x-y\, V^eK". 

Thanks to (H3), u is uniformly bounded by K. Moreover, (H4) and norm estimates of X in Lemma 2.1 ensure that 
u has the regularity properties given in the next Lemma. The proof is omitted since it is the same as the proof of 
Proposition 3.3 in [27], once we replace Lemma 3.1 of [27] by our Lemma 2.1. 

Lemma 2.2. Let us assume that g satisfies (H3) and (H4). Then there exists a constant L x > such that 

(2.6) |u(xi,i) — u(x 2 ,t)\ < L x \x\ — x 2 \, for any x\,x 2 GR and t £ [0, T]. 
Moreover, if the Levy measure satisfies (H2) ; then there exists a constant L t > such that 

(2.7) \u{x,tx) - u{x,t 2 )\ < L t |ti - i 2 | 1/2 , for anyt 1 ,t 2 6 [0,T] andxeR. 

The Lipschitz continuity of u(-,t) and the semi-Holder continuity of u{x, •) will be useful to show further regularity 
properties of u in the next three sections. 

For the optimal stopping problem, we define the continuation region C and the stopping region T> as usual: 

C = {(x,t) e R" x [0, T) : u{x, t) > g{x)} and T> = {(x, t) € R" x [0, T) : u(x, t) = g(x)} . 

2.2. The variational inequality. Intuitively, one can expect from the Ito's Lemma for Levy processes (see e.g. 
Proposition 8.18 in [9] pp. 279) that the value function u, defined in (1.2), satisfies the following variational 
inequality: 

min {(-d t - C + r) u(x, t), u(x, t) - g(x)} = 0, (x, t) € 1" x [0, T), 

(2.8) 

u(x,T) = g(x), 

in which the integro-diffcrcntial operator C, the infinitesimal generator of X, is defined via a bounded test function 
4> as 

(2.9) C4(x,t) ±C D <Hx,t)+I4>(x,t), with C D <t>{x,t) 4 a . j ^ t )—^- + J2b i (x,t)-^. 

i.j — l i—1 

Here A = (aij) nxn — \ cr(x, t)cr(x, t) T is a n x n matrix and the integral term 

. r T .Fid, 1 

v(dy) 



I<t>(x,t) 
(2.10) 



(x + y,t) - </)(x,t) - ^2y l Q—(x,t) l{| y |<i} 

i—l 



[<f>{x + y,t) - (j)(x, t)-y- V x (f>(x, t) l{\ y \<i}} v(dy). 



However, one does not know a priori whether the value function u is sufficiently regular (i.e., u e C 2,1 (R™ x [0, T))) 
to justify applying Ito's Lemma. Moreover, the integral term I(j)(x, t) is only well defined in classical sense when <p 
has certain regularity properties. It is sufficient to require </> to be a bounded function in C l (B t (x)), in which B t (x) 
is an open ball in R ra centered at x with some radius e g (0, 1), and that V x 0(-, t) to be Lipschitz in B e (x) uniformly 
in t, i.e., for t G [0,T) there exists a positive constant Lb such that \V x (j)(xi,t) — \7 x (j)(x2,t)\ < Lb\x\ — x 2 \, for 
xi,x 2 £ B e (x). Indeed, using these regularity properties of 4> we have that 

(2.11) I<f>(x, t) = I t <t>(x, t) + L c (f){x, t) < oo, where 
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(2.12) F<j>(x,t) = I [ ( j ) (x + y,t)-<j ) (x,t)]u(dy)-V x (f>(x,t)- yv{dy), 

\y\>t Jt<\v\<l 



(2.13) Ie<t>{x,t) = / {(j)(x + y,t)-cl>(x,t)-yV x 4>(x,t)]v(dy) 

J\y\<t 

= f f]y i (S lt i^(2i I t)-a B( ^ ) t))l/(*)< / L B \y\ 2 v{dy). 
J\y\<*l^{ J\y\<* 

(In (2.13), the second equality follows from the mean value theorem with some zt £ R" satisfying \zi — a; < |y|, while 
the inequality follows from the Lipschitz continuity of V x cf>(-, t). Note that J e< | J/ i< 1 \y\ v(dy) < i / e< ui<i 2/| 2j/ (^2/) < 
+oo, we obtain the inequality in (2.11).) 

Given the regularity of u in Lemma 2.2, it is not clear that u has the Lipschitz continuous first derivative to 
ensure that Iu is well defined in the classical sense. This will be addressed in Section 4. 

2.3. Viscosity solutions. We will introduce the viscosity solution concept following [27]. Let us define 

d(R n x [0, T}) 4 J £ C°(R" x [0, T]) : sup < +0O I . 

^ (i,i)eR»x[0,T] 1 + Fl J 

We adapt the notion of viscosity solutions used in Definition 2.1 of [27] into our context and give the following 
definition. (We assume that (H2) holds so that Icf>(x,t) is well defined for £ C 2ll (R" x [0,T]) n Ci(R" x [0,T]). 
Indeed, for 4> £ Ci(R™ x [0, T]), we have |</>(a; + y, <) — 0(cc, <)| < C(l + \y\) for some C independent of y. Therefore, 
f,, >e [<j>(x + y,t) - 4>{x,t)} v{dy) < +oo in (2.12) thanks to (H2).) 

Definition 2.1. (i) Any u £ C°(R™ x [0,T]) is a viscosity supersolution (subsolution) of (2.8) if 

(2.14) mm{-d t <f>- £(/) + ru,u(x,t) - g(x)} > (< 0), 

for any function (j) £ C 2a (R' 1 x [0,T]) n Ci(R™ x [0,T]) such that u(x,t) = (f>{x,t) and u{x,t) > (f>(x,t) (u(x,t) < 
4>{x, t) ) for all (x, t) £ R n x [0, T). 

(ii) u is a viscosity solution of (2.8) if it is both supersolution and subsolution. 

As in [27], we have the following result. 

Proposition 2.1. // the Levy measure v satisfies (H2), the value function u{x,t) is the unique viscosity solution 
of (2.8). 

2.4. Standing assumptions on the Levy Measure. In addition to (H2), throughout this paper, we assume 
that the Levy measure v has a density, which we denote by p(y), and there exists a positive constant M such that 

M 

(H5) p(y) < | a , for \y\ < 1 and some constant a £ [0,2). 

Remark 2.2. The Levy measures v , corresponding to Levy processes widely used in the financial modelling for the 
single asset case, satisfy (H5) with n = 1. 

In jump diffusions models where v is a probability measure, if the density p(y) is bounded, (H5) is satisfied 
with sufficiently large M. Examples of this case are Merton's model and Kou's model. On the other hand, if 
p(y) £ C°(B 1 (Q) \ {0}) and p(y) has a power singularity l/|j/| /3 with < /? < 1 at y — 0, (H5) is again fulfilled 
because \ y \i+ a > y^p- for any a > and \y\ < 1. 

Moreover, for Levy processes that are the Brownian motion subordinated by tempered stable subordinators, it 
follows from (4-25) in [9] that p has a power singularity l/|y| 1+2 ^ 7 with < j3 < 1, at y = 0. Therefore (H5) is 
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satisfied by choosing a = 2/3 and sufficiently large M. In particular, this class of Levy processes contains Variance 
Gamma and Normal Inverse Gaussian where (3 = or 1/2 respectively. 

Furthermore, for the generalized tempered stable processes (see Remark 4-1 in [9]) whose Levy density is p(y) = 

C xll C x 

\ y \ 1 + a - e 1{»<0} + y |i+Q + e ~ l{t/>o}; with a_,a+ < 2, (H5) is satisfied by choosing a = max{a_, a+, 0} 
and M = max{C_,C+}. In particular, CG MY processes in [8] are special examples of generalized tempered stable 
processes. In the similar manner, one can also check that the regular Levy processes of exponential type (RLPE) 
in [7] also satisfy (H5). 

3. Finite variation jumps and regularity in the continuation region 

In this section, we will analyze the regularity of u in the continuation region, when jumps of X have finite 
variation, i.e., 

(3.1) / (|y|Al)Kdtf)<+oo. 



It is clear that (3.1) holds when we assume (H5) is satisfied with < a < 1. In the main result of this section 
(Proposition 3.1), we show that u £ C 2,1 (C). 

Thanks to (3.1), the infinitesimal generator C can be rewritten such that its integral component has a reduced 
form, i.e., C(j)(x,t) = C f D (t)(x,t) + if cf>(x,t) in which C f D (f>{x,t) = £"i=i a^i^ + EiU h i - J|y|<i 2/ V ( d y) <W 
and 

(3.2) P'cb(t, x)= f [ct>{x +y,t)- <f>(x, t)} v{dy). 

Thanks to this reduced integral form and the Lipschitz continuity of u(-,t) (see Lemma 2.2), I^u(x,t) is well 
defined in the class sense. Indeed, it follows from (3.1) and (H2) that \l'u(x,t)\ < L \u(x + y,t) — u(x,t)\ v(dy) < 
L x J m \y\ v(dy) < +oo. Moreover, P 'u(x, t) is Holder continuous in its both variables. 

Lemma 3.1. Let £1 be any compact domain in R™ . Assume that the density p{y) satisfies (H5) with < a. < 1. 

(i) For any (xi,t), (x2,t) G O x [0,T], there exist constants Cn,/3 and Cq, independent of x\, xi and t, such 
that 

(3.3) when a = : \l^ u(x\,t) — I f u(x2,t)\ < C^,p\x\ — X2\ 1 ~ 13 , for any f3 £ (0,1); 

(3.4) when < a < 1 : \l f u(xi,t) - I f u(x 2l t)\ < Cn\xi 



X2\ 



(ii) For any (x,t\), (x,t2) G Q x [0,T] , there exist constants and Dq, independent of t\, t 2 and x, such 

that 



(3.5) when a = 0: \P u{x,ti) - I s u(x,t 2 )\ < Dn,f)\ti - t 2 \ 1 ^ , for any f3 £ (0,1); 

(3.6) when0<a<l: \l f u(x,h) - I f u(x,t 2 )\ < D n \h - t 2 | i ^- 

Proof. Our proof is motivated by Proposition 2.5 in [29]. See Appendix A for details. □ 

Now let us analyze (2.8) on a given compact domain inside the continuation region C. Let B be an open ball in 
l n such that B x (t 1 ,t 2 ) C C for some ti,t 2 £ [0, T). We will denote the closure of B by ~B. 

Instead of working with the nonlocal equation (2.8), note that V ' u is well defined in the classical sense, we 
consider the following local equation with the driving term I^u: 

(-d t - C f D +r)v(x,t) = I f u(x,t), for (x,t) £ B x [t u t 2 ), 
(3-7) _ 

v(x,t) = u(x,t), for (x,t) £ dB x [ii,t 2 ) U B x t 2 . 
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A viscosity solution of (3.7) is denned as follows. (See e.g. Definition 7.4 in [10], Definition 13.1 in [13].) 
Definition 3.1. Any v G C a (B x [iijta]) is a viscosity subsolution of (3.7) if 

(3.8) (-d t -£ f D +r)<t>{x,t) <I f u(x,t), for {x,t) 6 B X [t u t 2 ), 

(3.9) min j(-d t - C f D +r)4>(t,x) -I f u(x,t), v{x,t) -u{x,t)} < 0, for (x,t) G dB x [t 1 ,t 2 )UB x t 2 

for any function <fi G C 2,1 (K™ x [t\, t 2 }) such that <fi(x, t) = v(x, t) and </>(i, t) > v(x, t) for any (x, t) G M™ x [ti, t 2 ]. 
The supersolution is defined analogously. As usual, v is a viscosity of (3.7) if it is both a subsolution and a 
supersolution. 

Based on Proposition 2.1 and an equivalence between Definition 3.1 and another notion of viscosity solutions 
for a nonlocal equation, we have the following result, whose proof is listed in Appendix A. 

Lemma 3.2. If the Levy measure v satisfies (H2) and (H5) with < a < I, then u is a viscosity solution of (3.7). 
Now let us assume that the diffusion component of X is nondegenerate i.e., there exists A > such that 

n 

(H6) aij{x,t)C£ j > A|£| 2 , Va;,£ £ R n ,t > 0. 

Additionally, we assume 

(H7) aij(x,t),bi(x,t) and r(x,t) are continuously diffcrentiable in both variables on M™ x [0, T] for i,j < n. 
We are ready to show the regularity of u inside the continuation region. 

Proposition 3.1. Assume that the Levy measure v satisfies (H2) and (H5) with < a < 1, moreover (H6) and 
(H7) are satisfied. Then u £ C 21 (C). 

Proof. First it follows from Lemma 2.2 that the boundary and terminal values of (3.7) are continuous on dB x 
[ti,t 2 ) U B x t 2 . Moreover, the driving term Pu{x,t) is uniformly Holder continuous in B x [ti,i2] thanks to 
Lemma 3.1. Therefore, (3.7), whose coefficients satisfying (H6) and (H7), has a unique classical solution u*(x,t) G 
C 2,1 (B x (ti,t 2 )) (see Theorem 9 in [14] pp. 69). Since u* is already a classical solution, u* is also a viscosity 
solution of (3.7). Now, it follows from the comparison theorem for viscosity solutions (see e.g. Theorem 7.5 in 
[10]) that u(x,t) = u*(x,t) for (x,t) G B x (ti,t 2 ). Therefore u G C 2 ^(B x (t\,t 2 )). The statement follows since 
B x (ti,t2) is chosen arbitrarily in C. □ 



4. Infinite variation jumps and the regularity on the whole domain 

In this section, we shall study the regularity of the value function on the whole domain. Moreover, jumps in 
this section may have infinite variation (i.e., (3.1) may not be satisfied). We give the main result of the paper in 
Theorem 4.1, which states that u G Wp' 1 (B), for any compact domain B C M. n x (0,T) and p > 1. The proof of 
this result is given in Section 5. There are two important corollaries to Theorem 4.1: In Corollary 4.1, we show 
that the smooth fit property holds; in Proposition 4.1 we show that u G C 2,1 (C). We start by developing some 
properties of the integral operator / in Lemma 4.1. These properties will be crucial in our proofs. 
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4.1. Notation. 

Definition 4.1. Let fi fee a domain in R™ ; Qt = fi x (0, T) and Qt be the closure of Qt- 

C (Qt) denotes the class of continuous functions v on Qt with continuous classical derivatives dtV, d x iV and 

For any positive integer p > 1, Wp' (Qt) is the Banach space consisting of the elements of v € L p (Qt) having 
generalized derivatives of the form dtv, d x iV and d x i x jV for i,j < n and a finite norm \\v\\ w ^,i ^Q T y (Please refer 
to pp. 5 of [22] for the definition of \\ ■ \\ w 2.i .) 

In particular, W^^JJQt) is the Banach space consisting of functions whose W^-norm is finite on any compact 
subsets of Qt- 

For any positive nonintegral real number a, H a ' a ^ 2 (Qt) is the Banach space of functions v that are continuous in 
Qt, and have continuous bounded classical derivatives of the form d^d x v for2r + s < a, which satisfy < oo. 

(Please refer to pp. 7 of [22] for the definition of \\ ■ \\' a ' ■) 

On the other hand, H a (fi) is the Banach space whose elements are continuous functions v(x) on fi that have 
continuous bounded derivatives up to order [a] and satisfy \\v\\~ < oo. 

4.2. The integral term. When the jumps of X have infinite variation, i.e., (3.1) is not satisfied, the integral 
term cannot be reduced to the form in (3.2). Therefore, throughout this section we need to work with the integro- 
diffcrcntial operator L and its integral part J in the form of (2.9) and (2.10). However, given the regularity 
properties of u in Lemma 2.2, it is not clear that Iu is well defined in the classical sense. (See the discusion after 
(2.10).) Nevertheless, we shall show in the following lemma that given sufficient regularity properties for the test 
function <f>, I<fi(x, t) is Holder continuous in both variables. Later in this section, we will prove that the value 
function u does have these regularity properties to guarantee Iu well defined in the classical sense. 

Let fi be a compact domain in R", fi 5 = {x E R n : x E Bg(y) for some y € fi} for some S > 0. For s G (0,T], 
let us denote Q7 = fi x [0, s] and <?J = TF x [0, s}. Moreover, we denote D s = R" x [0, s\. 

Lemma 4.1. Let us assume that the Levy measure satisfies (H2) and (H5) with a G [1,2). 

(i) Assume that <f> satisfies max R n x [ 0jS ] \<f>\ < oo, max K n x [ . s ] |V X <^| < oo, and that there exists L t > such 
that \4>(x,h) - (f>(x,h)\ < Lt \h - i 2 1 1/2 for any x G R and ti,t 2 G [0,s]. If $ E H 13 '? (pfj for some 

f3 E (a, 2), then Iu E H^~ , ^~ (Q s ). Additionally, there exists a constant Cn > 0, depending on fi, a, j3 
and T , such that 

(4.1) < C n ( max |0| + max |V^| +L t + \\</>\\W 

(ii) If 4> E H l3 '^(D s ) for some f3 E (a, 2), then 1(f) E '^~z~(D s ). Moreover, there exists a constant C, 
depending on a, (3 and T, such that 

(4.2) imi^ <CU\\ { £. 

Proof. The proof is similar to the proof of Lemma 3.1. See Appendix A. □ 

Remark 4.1. When the Levy measure v is a finite measure on R n , the integral form J Rn (f>(x + y, t) v(y) has the 
same regularity as <p(x,t); see [31]. When the Levy measure has a singularity, as we have seen in Lemma ^.1, the 
mode of continuity of I <f> decreases compared to the mode of continuity of <fi. Moreover, as we have seen in (4.1), 
the Holder norm of I <fi depends on the Holder norm of <j> on a slightly larger domain. This extension of domains 
will introduce a technical difficulty in estimating the Sobolev norm of u. This estimation will be carried out in the 
following section. 
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4.3. Solutions in the Sobolev sense. In this subsection, we shall give the main result of this paper. Compared 
to Section 3, we need some stronger assumptions on coefficients. Instead of (H7), we assume that 

(H7') a ij (x,t),b(x,t),r{x,t) e H e ^(W l x [0,T]), W G (0, 1) and i, j < n, and r(x, t) > 0. 

Moreover, there exist positive constants A such that 

n 

(H6') A|£| 2 < Yj Oijfat)?? < A|£| 2 , W(x,t) el"x [0,T] and £ G R n . 

In addition to (H3) and (H4), we assume that there exists a positive constant J such that 

d 2 

(H8) tt^9>-J, mS'fW 1 ), for any direction ry G R™, 

in which d/drj is the directional derivative, and the inequality is understood in the distributional sense. 

Let £ e be the standard mollifer (see [7] pp. 629 for its definition and properties). Consider the mollified sequence 
9 £ — .9 * C c f° r e G (0, eo)- Here eo is a positive constant less than 1. First, it follows from (H8) that 

n 

(4.3) J2d 2 xixj g*( X )et; j >-J\e, V^zgR". 
It is clear that 

(4.4) each g e (x) G # 2+ ^(R n ) W G (0, 1). 

Additionally, (H3) and (H4) imply that there exist positive constants K and L independent of e such that for all 

x G R™ 

(4.5) < g e (x) < K, 



(4.6) |V 5 e (x)| < L. 

Now we are ready to state main result of this paper. 

Theorem 4.1. // (H3), (H4), (H6), (H7 1 ), and (H8) are satisfied, and the Levy measure v satisfies (H2) and (H5) 
with a G [0, 2), then u G W 2,1 {B p {xq) x (0,T — s)) /or am/ integer p G (1, 00), p, ,s G M+ anrf x'o G R™. Moreover, 
u solves (2.8) /or almost every point in R" x [0, T]. 

Before wc prove this result in Section 5, let us list some of its corollaries. 

Corollary 4.1. // the assumptions in Theorem J^.l are satisfied, then for any p, s > and xq G R™ 

(i) u G HP<%(B p (x ) x [0,T - s]) where /3 = 2 - ^±2 > q In particular, V x u G C(R" x [0,T)). Therefore, 
the smooth-fit property holds. 

(ii) // the Levy measure v satisfies (H5) with a G [1,2), then Iu is well defined in the classical sense in 
B p (xo) x [0,T). Moreover, Lu G H^'^ (B p (x ) x [0,T - s]) /or some /? G (a, 2). 

Proof, (i) Combining Theorem 4.1 and the Sobolev Inequality (see e.g. Lemma 3.3 in [22] pp. 80), we have 
u G (B p (xq) x [0,T— s]), where /? = 2 — > 0. Choosing sufficiently largep such that (3 > 1, the continuity 
of Vxii follows from Definition 4.1. 

(ii) Let us choose p sufficiently large so that f3 > a. Now, the proof follows from (i) and Lemma 4.1. □ 
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Thanks to Corollary 4.1 (ii), we can consider the following boundary value problem with the driving term I u: 

{-d t - C D +r)v(x,t) = Iu(x,t), for (x,t) £ B x [ti,t 2 ), 
(4.7) _ 

v(x, t) = u(x, t), for (x, t) £ dB x [ti,t 2 ) U B x t 2 , 

where Bx (ti,t 2 ) C C is the bounded domain as in (3.7). The viscosity solution of (4.7) is defined as in Definition 3.1, 
after replacing C^ D and V by Co and I, respectively. 

The following relation between the solutions in the Sobolev sense and the viscosity sense shows that u is a 
viscosity solution of (4.7). See Corollary 3 in [25] or Theorem 9.15 (ii) in [20] for its proof. 

Lemma 4.2. If u £ Wp' 1 (B x (ti,t 2 )) for p > n + l satisfies (4.7) at almost every point in B x (ti,t 2 ), then u is 
the viscosity solution of (4.7) in the sense of Definition 3.1. 

Thanks to Corollary 4.1, Lemmas 4.1 and 4.2, the arguments in the proof of Proposition 3.1 now works for the 
infinite variation jump case. 

Proposition 4.1. If the assumptions of Theorem 4-1 are satisfied, then u £ C 2 l {C). 

Proof. Corollary 4.1 (ii) tells us that Iu(x,t) is Holder continuous in both its variable. Moreover, u is a viscos- 
ity solution of (4.7) thanks to Lemma 4.2. The rest proof follows the same line of arguments in the proof of 
Proposition 3.1. □ 

In our future work, we will investigate the regularity of the free boundary curves and extend our results in [5]. 

5. Proof of Theorem 4.1 

Because the jump may have infinite variation, the proof of Theorem 4.1 needs to conquer several technical 
difficulties. We will carry the proof of Theorem 4.1 in a series of lemmas and point out the difficulties along the 
way. 

Let us first define v(x, t) — u(x, T — t) for (x, t) £ R n x [0, T]. It is natural to expect that v solves the following 
variational inequality 

min{(9 t - C D - I + r) v(x, t), v(x, t) - g(x)} = 0, (x,t) £ R x (0,T], 

(5.1) 

v(x,0) = g(x). 

We will establish Theorem 4.1 by using the penalty method, which constructs a sequence of approximating functions 
each of which solves (5.2). First, in Lemma 5.2, we find a nice enough solution, v t , to each penalty problem. In 
Corollary 5.1 we give a uniqueness result for each of these penalty problems. Second, we analyze the properties of 
the value functions of the penalty problems in Lemmas 5.4, 5.5, 5.6, and Corollary 5.2. These are used to show 
that the W^-norm of v e is bounded uniformly in e in Corollary 5.3. In order to establish the latter result, we also 
prove a -norm estimate for the solutions of parabolic integro-diffcrcntial equations in Theorem 5.1. We show 
in Corollary 5.4 that the weak limit of {v e }, which we denote by v*, solves the variational inequality and has a 
finite Wp' 1 norm. This result along with Proposition 5.1 concludes the proof of Theorem 4.1. 

In the following, we will only carry out the proof of Theorem 4.1 for the infinite variation jump case, i.e., the 
Levy measure v satisfies (H5) with 1 < a < 2. Since the integral operator has the reduced form I* in (3.2) for the 
finite variation jumps, the proof of < a < 1 case in Theorem 4.1 will be similar and easier. 

Motivated by Lemma 3.1 in [15] pp. 24 and [31], we will study the following penalty problem for each e £ (0, eo): 

(dt-£ D -I + r)v € (x,t)+p e (v £ -5 £ ) = 0, (x,t) £ R" x (0,T], 

(5.2) 

« £ (x,0) =/(*), 
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in which {.g e } e6( - eo ) is given after (H8). Here, the penalty term p e (y) £ C°°(R) is chosen to satisfy following 

properties: 

(5.3) 

{i)Pe(y)<0, (ii)p e (y)=0i£y>e, (in)p e (0) = -nAJ- - \r\<®K- J [ \y\ 2 u{dy) - K [ v(dy), 

J\y\<l J\y\>i 

(iv)p' e (y)>0, (v)p e (y)<0 and (vi) ]imp e (y) = \ °' V > ° . 

e|0 — oo, y < 

The constants A, K, L and J come from (H6), (4.5), (4.6), and (H8), respectively. Additionally, =max K ™ x [ 0; T] \b( 
and \r\(°' = max Rnx [ .T] |r(x, i)| a re finite due to (H7'). Moreover, p e (0) is also finite thanks to (2.3). It is also 
worth pointing out that p c (0) is independent of e. These properties of p t will be useful in the development of our 
next few results. In particular, (5.3) (hi) is essential for proofs of Lemma 5.6 and Corollary 5.2. 
Let us recall the Schauder fixed point theorem (see e.g. Theorem 2 in [14] pp. 189). 

Lemma 5.1. Let © be a closed convex subset of a Banach space and let T be a continuous operator on such 
that TO is contained in © and TO is precompact. Then T has a fixed point in 0. 

For each e £ (0, eo), we will show that the penalty problem (5.2) has a classical solution via the Schauder fixed 
point theorem. Let us recall D s = M" x [0, s]. 

Lemma 5.2. If the Levy measure v satisfies (H2) and (H5) with 1 < a < 2, then for any e £ (0, eo) and fj G (a, 2), 
(5.2) has a solution v e 6 H 2+i ^- 1+! ^ £L (D T ). 

Proof. We will first prove that (5.2) has a solution on a sufficiently small time interval t £ [0, s] via the Schauder 
fixed point theorem. Then we will extend this solution to the interval [0,T]. 

Let us consider the set 6 = (v £ H^'^(D S ) with its Holder norm < Uoj, where positive constants s and 

Uq will be determined later. It is clear that is a bounded, closed and convex set in the Banach space H^'i (D s ). 
For any v £ 0, consider the following Cauchy problem for u — g t : 

(d t -C D +r){u~ g e )(x, t) = Iv(x, t) -p e (v- g e ){x, t) + (C D - r) g e (x), (x, i)elx (0, s], 

(5.4) 

u{x,0) -g e (x) = 0. 

Via the solution u of (5.4), the operator T can be defined as u = Tv. Let us check the conditions for the Schauder 
fixed point theorem are satisfied: 

1. Tv is well defined. Since v £ H@'^(D S ) and (3 £ (a, 2), it follows from Lemma 4.1 (ii) that Iv £ 
H^^{D S ) with 

(5.5) II-^IId 2 < C \\v\\ jj , for some constant C > independent of s. 

On the other hand, we can check that p e (v — g t ) £ H^s~ (D s ). Indeed, p e (v — g e ) is bounded in D s , since both 
v 1 g t £ HP'2(D S ) (see (4.4)) and p e (y) £ C°(R). Additionally, for any x u x 2 £ K", t £ [0, s] 

|p e (« - 5 e )(a;i,i) -p e (w - 9 £ ){x2,t)\ < max |p e (w - g e )\ \(v - g £ )(x 1 ,t) - (v - g 6 ){x 2 ,t)\ < C\xi - x 2 \- 

Here max£> 5 \p e {v — g e )\ is finite, which also follows from the boundness of v — g e and p e £ C 1 (R). The positive 
constant C depends on maxc s \p e (v — g e )\ and the Holder norms of v and g e . The Holder continuity of p e (v — g e ) 
in t can be checked similarly. Furthermore, {Co — r)g e (x) £ (D s ) as a result of (4.4). Therefore, 

thanks to (H6) and (H7'), it follows from Theorem 5.1 in [22] pp. 320 that (5.4) has a unique solution u — g e £ 
H 2+ ^r> 1+ ^(D s ). Note that g e £ H 2+S ^' 1+e ^ (D s ) (see (4.4)). As a result u = Tv £ H 2+M ^ i ' 1+i ^ 1 (D s ). 
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2. T0 C 6. For u = Tv, appealing to Lemma 2 in [14] pp. 193, we obtain that there exists a positive constant 



Ap, depending on /?, such that 
(5.6) 



\Iv\r + \\ Pt (v-f)\\^ + \\{£ D -r)f\\ 



(0) 



<A p Csi\\v\\f +A, 



where 7 = ^-5^, C is the constant in (5.5) and A is a sufficiently large constant dependent on Hc^Hr,^ for some 
i G (0,1). Let s be such that r = A, 3 Csi < 1/2 and let U = max{ T ^r, 2 ||ff e ||^}. Note that \\v\\ { £] < U . Now 
it follows from (5.6) that 



(5.7) 



NlK<ll«-0 e ||g? 



W\\^<rU a + A 



U ° <r tt _i_ !- 2t tt j. U ° tt 
-j <tU + — — U + y = u o- 



Therefore, u = Tv G 0. 

3. T0 is a precompact subset of H^'t(D s ). For any 77 € (/3, 2), similar estimate as (5.6) shows that for 
any v G O, we have ||Tu||^ < U± for some constant U± depending on Uq and s. On the other hand, argument 
similar to Theorem 1 in [14] pp.188 shows that bounded subsets of H^-i (D s ) are precompact subsets of H^'i (D s ). 
Therefore, T0 is a precompact subset in H^'^{D S ). 



4. T is a continuous operator. Let v n be a sequence in O such that lim 



I 03) 



0, we will show 



lim 



3 \\Tv n — Tv\\p = 0. From (5.4), w = Tv n — Tv satisfies the Cauchy problem 

(d t -C D +r) w(x, t) = I(v n - v)(x, t) - [ Pe {v n - <? e ) - p e (v - g e )} , {x, t) G R" x (0, s] 
w{x,Q) = 0. 

It follows again from Lemma 2 in [14] pp. 193 that 



\Tv n - Tv 



i(/3) 
D„ 



< A 



/3S' 



|/K-«)|| (0) + ||p £ K-.g e )-p e («-ff e 



C||u„ - + max p e (v n - g e ) \\v n - v\\ 



D s .n 



|(0) 



|(0) 



as n-> 00. 



As a result of Steps 2. - 4. and the Schaudcr fixed point theorem, we obtain a fixed point of the operator T 
in H I3 ^(D S ). We denote this fixed point by v e . Moreover, it follows from the result in Step 1 that w e = Tv e G 
H 2 +^++^{D s ). 

Finally, let us extend v e to the interval [0,T]. Choosing any p G (0,T — s), we replace g e (-) by v e (-,p) in 



is finite, we 



(5.4). Note that the choice of s in 2. only depends on j3 and C, but not on p. If |u e (-, /9)||m»» 

(21 ) 

can choose a sufficiently large Uo, depending on ||v e (-, p)\\ R n 2 , such that (5.7) holds on [p, p + s], moreover 
||w £ (-, /O + s)||^ 2 „ + 2 ' is finite thanks to the result after 4.. Noticing that ||g e ||R n + ^ is finite for any I G (0, 1), one 
can extend the time interval by s each time, until the time interval contains [0, T]. Therefore we have the statement 
of the lemma. □ 

Thanks to the definition of the Holder spaces, Lemma 5.2 also tells us that v e is bounded in Dt- In order to 
show that v e is the unique bounded classical solution of the penalty problem (5.2), we need the following maximum 
principle for the parabolic integro-differential operator. The proof of it is provided in Appendix A. (See Lemma 2.1 
of [31] for a similar maximum principle, where v is assumed to be a finite measure on R.) 

Lemma 5.3. Let us assume that a.ij(x,t), bi(x,t) and c(x,t) are bounded in R™ x [0,T] with A = (aij) nxn 
satisfying Y17j=i a ij{ x i *) > f or an y £ € E n \ {0}, moreover c(x, t) > and the Levy measure satisfies (H2). 
Lfv G C°([0,T] x M")nC 24 ((0,T] x R") satisfies (d t - Cv — I + c(x, t)) v(x, t) > in Rx (0,T] and there exists a 
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sufficiently large positive constant m such that v(x,t) > —m for (x,t) E R™ x [0,T]. Then v{x, 0) > implies that 
v{x, t) > for (x, t) E R" x [0, T] . 

As a corollary of this maximum principle, the bounded classical solution of the penalty problem (5.2) is unique. 

Corollary 5.1. For each e E (0, eo); the penalty problem (5.2) has a unique bounded classical solution. 

Proof. Let us assume v\ and V2 are two bounded solutions of (5.2). Then V\ — v 2 satisfies 

(d t -C D -I + r) (vi - v 2 ) + Pe (vi - g £ ) - Pe (v 2 - g e ) = 0, (x, t) E R n x (0, T] , 

(5.8) 

(«i-«2)(a?,0) = 

On the other hand, it follows from the mean value theorem that p e (v\ — g e ) — p e {v 2 — # e ) = P € {y){vi ~ V2) for 
some y E R™. Moreover, p t (y) is bounded, say by M, thanks to the fact that p e E C 1 (R) and V\, v 2 and g e are all 
bounded. Now applying Lemma 5.3 to the equation (5.8) and choosing c = r + M > (see (5.3) (iv)), we have 
vi(x,t) > v 2 (x,t) for (x,t) E M" x (0, T]. The other direction of the inequality follows from applying the same 
argument to v 2 — Vi . □ 

Applying Lemma 5.3, we will analyze some universal properties of v e for all e E (0, eo) in the following three 
lemmas. 

Lemma 5.4. 

< v e (x,t) <K+1, for (x,t) el"x [0,T]. 
Proof. Since the proof is similar to the proof of Lemma 2.2 in [31], we give it in the Appendix A. □ 
Lemma 5.5. 

\d x kV € (x,t)\ < L, for (x,t) E R n x [0,T],1 < k < n. 
Proof. The proof is similar to the proof of Lemma 2.4 in [31]. □ 
Lemma 5.6. For any e E (0,e ) ; v e (x,t) > g e (x) on R n x [0,T]. 
Proof. Let us first show that Ig e (x) is uniformly bounded from below. Indeed, 





wl<i 



g e (x + y)- g e (x) - ^ y *— g *(x) 

i—1 



v{dy)+ I [g e (x + y)-g e (x)]u{dy) 
|y|>i 



/ v(dy) / dz{l-z) V yY——g*( x + zy)+ \ 
4l<i Jo flti dxdxJ Av\>i 



[g t {x + y)-g'{x)]v{dy) 



(5.9) Jm - " u «=i 

> 



u{dy) / dz(l - z) (-J\y\ 2 ) - K / u{dy) 

y\<l JO J\V\>1 



>-J \y\ 2 v(dy)-K V {dy), 

J\y\<l J\v\>l 

where the first inequality follows from (H8) and (4.5). 

On the other hand, thanks to (H6) and (H8), Y^ij a ij( x jt)d 2 i x jg e (x) is also bounded from below. Note that 
<Hj(x,t) d% ixj g e (x) = tr(AH(g e )), where H(g e ) is the Hessian of g e , i.e., H(g e ) i3 = d* ixj g e (x). It follows from 
the first inequality in (H6) that A is a positive definite matrix. Then there exists a nonsingular matrix C such that 
A = CC' . Therefore tr(AH(g e )) = tr{CC' H(g € )) = tr(C' H(g e )C). Moreover, (H8) and (H6) give us that 

{CO H{g e ) {CO > -J U'C'CA = -J U'At) > -JAICI 2 , E R". 
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Hence C H(g e )C + JAI n is a non-negative definite matrix. As a result, we have tr (C H{g' i )C) + nJA 
tr (c'H(g e ) C + JAI n ) > 0, which implies 



(5.10) MM) d^dJ 9 ^ x) = tr ( AH{gC)) - ~ nJA - 

Thanks to (5.9) and (5.10), we can bound (dt — £d — I + r) g e {x) from above. Indeed, 
(d t -£ D -I + r)g e (x) 

" d 2 " d 

= - MM) dx i dxj 9 £ ( x ) -^h{x,t) — i g t (x)+r(x,t)g t (x) - Lg c { 3 

(5.11) l ' J=1 1=1 

< nJA + \b\( Q) L + \r\( 0) K + j[ \y\ 2 v{dy) + K [ v{dy) 

J\y\<l J\V\>1 

= ~Pe(0), 

where the second equality follows from (5.3) (hi). 

Now we will show that v e > g e using Lemma 5.3. It follows from (5.11) that 

(d t -C D -I + r) (u e - <f ) = -p e (v* -g*)-(d t ~£ D -I + r)g ( 

> -p e (v e -g e )+Pe(0). 

The last equation together with the mean value theorem implies that 

(5.12) (d f -£ D -I + r + p'M) (v e - 5 e ) > 0, 



for some y € R™. Therefore the statement of the lemma follows applying Lemma 5.3 to (5.12) and choosing 
c = r + p' e (y) > 0. □ 

As an easy corollary, the penalty terms are uniformly bounded. 

Corollary 5.2. p e (v e — g e ) is bounded uniformly in e £ (0, eo). 

Proof. Thanks to Lemma 5.6 and (5.3) (i) and (iv), we havep e (0) < p e (v e — g e ) < 0. The statement follows noticing 
that p e (0) (in (5.3) (hi)) is independent of e. □ 

Thanks to Lemmas 5.2, 5.4, 5.5 and Corollary 5.2, we can apply the following W^-norm estimate for the 
parabolic integro-differential equation to each solution v t of the penalty problem. 

Since the proof of the following theorem is technical and independent of the penalty problem, we will perform it 
in Section 6. 

Theorem 5.1. Let us assume the Levy measure satisfies (H5) with a € [0,2), if v is a W p ' loc solution of the 
following Cauchy problem for some positive integer p, 

(dt-C D -I + r)v = f(x,t), (x,t) 6 1" x (0,7], 

(5.13) 

v(x,0) = g(x), 

where the coefficients satisfy (H6), (H7') and f £ L p j oc (M. n x (0,T)), moreover \v\ is bounded on R n x [0, T] and 
\V x v\ is bounded on any compact domain o/R™ x [0, T] . Then for any domain B p (xq) x (s, T) with p > 0, s e (0, T) 
and x 6 R™ 



(5.14) \H wi , HBp{xo)x{StT)) <C s 
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for some positive constant C$ and S < s. 

Remark 5.1. The existence of the W^ 1 solution for (5.13) was ensured by Theorem 3.2 in [6] pp. 234- However, the 
norm estimation was not given there. On the other hand, since the integral operator I is non-local, it is important 
to study the Cauchy problem (5.13) on the entire domain K" x [0, T] . Otherwise, for the Cauchy problem on bounded 
domains ofW 1 x [0,T] with some boundary conditions, W?' solutions are not expected in general, see [17] for a 
counterexample. 

A Wp' -norm estimate, similar to (5.14), for the parabolic integro- differential equation was proved in Theorem 3.5 
in [1C] pp. 91. However, the estimation in [16] requires the jump restricted in a bounded domain, i.e., if x G Q 
where £1 is a bounded domain in W 1 , the jump size z(x), which is state dependent, can only be chosen such that 
x + z(x) G fl (see (1.54) in [16] pp. 63). However, this restriction is not satisfied in our case, where the jump size 
is unbounded and independent of the state variable x. 

Applying Theorem 5.1 to each penalty problem (5.2), thanks to Lemmas 5.2, 5.4, 5.5 and Corollary 5.2, we have 
the following corollary. 

Corollary 5.3. // the assumption of Theorem 4-1 are satisfied, then for any domain B p (xq) x (s,T) with p > 0, 
s G (0,T) and xq G M", ||w e || w/ 2,ij B (a; )x(sT)) are bounded uniformly in e G (0,eo) for any integer p G (l,oo), i.e., 
there is a constant C independent of e such that 

( 5 - 15 ) \\ ve \\w 2 p -\B p ( X0 )x(s,T)) < C - 

Proof. It follows from Lemma 5.2 that v e G Wpf oc {R n x (0, T)). Thanks to Lemmas 5.4 and 5.5, both max Rn x [o,t] \ ve \ 
and maxfnxjQ^] |V x u e | are bounded uniformly in e. Moreover, it follows from Corollary 5.2 that / = — p e (v e — g f ) 
is also bounded uniformly in e. Therefore, (5.15) follows from (5.14). □ 

Remark 5.2. Theorem 5.1 is essential for the proof of Corollary 5.3. However, having infinite variation jumps 
presents two technical difficulties to the proof of Theorem 5.1. First, as we shall see in Lemma 6.1, once the Levy 
measure has a singularity, the L p -norm of Iv e depends on the W^ 1 -norm of v e . Therefore, one could not consider 
Iv e as a driving term directly and use the classical W^ 1 -norm estimate for parabolic differential equations (without 
the integral term) to bound the Wp' 1 -norm ofv e by the L p -norm of I v e . On the other hand, when the Levy measure 
is a finite measure as in [31], the L p -norm of Lv e only depends on the L°° -norm of if. Therefore, Lemma 2.6 in 
[31] follows from the classical Wp ' -norm estimate for parabolic differential equations, i.e., the W p ,l -norm of if is 
bounded by the L°° -norm of v e . 

Second, as we have seen in Remark 4-1 and we shall see it again in Lemma 6.1, the regularity of Iv e actually 
depends on regularity of v e on a larger domain. This extension of the domain is another technical difficulty we face 
in the proof of Theorem 5.1, because the extension of domains implies that W^' 1 -norm of v e on a bounded domains 
depends on its W^-norm on a slightly larger domain. 

To conclude this section, in the following theorem we will find a limit v* of the sequence {u c } £6 (o ;£o ) and show 
that it is the value function v defined at the beginning of this section. 

Corollary 5.4. Let us assume that the assumptions we made in Theorem 4-1 are satisfied. Then for any s,p>0 
and Xq G R™, there exists a subsequence {e/c}fc>o such that v ek converges uniformly to the limit v* uniformly in 
B p (xq) x [s,T] as — } 0. Moreover, v* solves the variational inequality (5.1) for almost every point in R™ x [0,T] 
and v* G Wp' 1 (B p (xo) x (s,T)) for any integer p G (1, co). 
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Proof. Combining Corollary 5.3 and the fact that Wp' 1 is weakly compact, we can find a subsequence {e k } with 
e k -> and a function v* 6 Wp' 1 (B p (xo) x such that 

(5.16) ^ -a v* in W^Bpfco) x (s,T)). 

Here " — " represents weak convergence. Refer to Appendix D.4. in [12] pp. 639 for its definition and properties. 
The rest of the proof is the same as proof of Theorem 3.2 in [31]. It confirms that v* solves the variational inequality 
(5.1) for almost every point in W 1 x [0, T]. □ 

Next, we present a verification theorem relying on the stability and uniqueness of viscosity solutions adapting 
the arguments in the proof of Theorem 1 in [3] . 

Proposition 5.1. The function u*(x,t) = v*(x,T — t) is equal to the value function u. 

Proof. We will show that u* is a viscosity solution of (2.8). Then the fact that u* = u is a consequence of 
Proposition 2.1. 

First, let us show that u* is a viscosity subsolution of (2.8). Consider a point (x,t) and a test function <j> e 
C 2ll (R™ x [0,T]) n d(R n x [0,T]) such that u*(x,t) < 4>{x,t) = u*(x,t) for any (x,t) el"x [0,T). 

Let us assume that u*(x,t) > g(x), otherwise there is nothing to prove. Therefore, we need to prove that 

{d t -£ D -I + r)4>{x,t) <0, 

for such a point (x, t). In fact, it is enough to prove that for any test function 4> 6 C 2,1 (IR™ x [0, T}) nCi(R™ x [0, T]) 
such that u*(x,t) < <fi(x,t) = u*(x,t) for any (x,t) E B(x,t;S). 

(5.17) (dt4>-C D <t>-h4>-I 5 [u*,<t>}+r<l>) (x,t) < 0, 
in which Is is defined in (2.13) and 

l s {u*,0}(x,t)= [ [u*(x + y,t)-u*(x,t)}v(dy)-V^(x,t)- f yv(dy); 

J\y\>6 JS<\y\<l 

see e.g. Proposition 1 of [3]. 

Without loss of generality we can assume that (x,t) is the strict maximum of the function u* — <fi on B(£,cc;<5), 
otherwise the test function can be appropriately modified. Then, thanks to the uniform convergence of u tk (x, t) = 
v tk (x, T — t) to u* , and the continuity of the functions in this sequence, we have that for large enough k the function 
v tk — 4> attains its strict maximum over B(x,t;5) at (xk,tk) € B(x,t;S). Moreover, (xk,tk) — > (x,t). Since v £k is a 
classical solution of (5.2) (see Corollary 5.1), it is also a viscosity solution. Hence (again using something similar 
to Proposition 1 of [3]), 

(5.18) {d t <j> ~ L D <f> - l s( j> - I 5 [u e \<j>] + r<£) (x kl t k )+ Pek {u e -{x k ) - g ek (x k j) < 0. 

First, since u*(x, t) > g(x) and that u tk {x k ) — g tk {x k ) — > u*(x,t) — g(x) the penalty term in (5.18) goes to 
0. Second, using (2.13) and the dominated convergence theorem we have that Is4>(x k ,t k ) —> Is<f>(x,t). Third, 
using the fact that u tk is uniformly bounded (see Lemma 5.4) we can exchange the limit and integration in 
lim.fc_nx3 I s [u tk , (j>](xk, t k ). Furthermore, using the uniform convergence of u ek to u* on compacts, we conclude that 
I s [u ek , <p](x k ,t k ) — > I s [u*, <ft](x, t). Therefore, passing to the limit in (5.18) we obtain (5.17). 

Second, we will prove that u* is a viscosity supersolution of (2.8). Since u ek > g e by Lemma 5.6, it follows that 
u* > g. Therefore, we only need to show that 

{d t -C D -I + r)4>{x,t)>Q, 
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for a test function <f> g C 2 ^(R n x [0,T])nCi(R" x [0,T]) such that u*(x,i) > </>(x,t) = w*(ar,i) for any (x,t) 6 
]R n x [0,T). The rest of the proof of the supersolution property follows the line of arguments we used for the 
subsolution property. 

□ 



Theorem 4.1 follows from Corollary 5.4 and Proposition 5.1. 



6. Proof of Theorem 5.1 

For notational simplicity, the constant C denotes a generic constant in different places. Moreover, the center xq 
of the ball B p {xq) will not be noted in the sequel. For any positive integer p, let us first estimate the i p -norm of 
the integral term Iv. 



Lemma 6.1. If the assumptions of Theorem 5.1 are satisfied, then for any rj > ; there exists a positive constant 

C such that 

(6.1) 

' (1 + ?/ 1 -"), a^l 



\\Iv\\l p (B p (x )x(s,T)) <Cv a \\v\\ w 2,i iBp+v{xo)x{sT) y 



C I max Id + max |Vx«| )• 

>x[s,T] B p+1 (x )x[s,T] 



(1 - log rf), a = l 



Proof. Let us break the integral into three parts. 



\Iv(x,t)\ 



[v(x + y,t)- v{x, t)-y V x v(x, t)l{\ y \<i}] v{dy) 



< 



. .1 n 

/ v{dy) / dz(l -z)J2 
J\v\<v Jo f , =1 



2/V 



d 2 



dx l dxi 



v(x + zy, t) 



v(dy)\v(x + y,t) -v(x,t) -y V x v(x,t)\ + / v(dy)\v(x + y,t) - v(x,t)\ 
1<\v\<l J\v\>l 



< V / \y\ 2 v{dy) / dz 

ij=lJ\v\<r> J 



if- 



dx t dx J 



v(x + zy,t) 



u(dy) \v(x + y,t) -v(x,t) -y ■ W x v(x,t)\ + / u(dy)\v(x + y,t) - v(x,t)\ 
v<\v\<i J\y\>i 

71 

: I^(x,t)+I 2 {x,t)+I 3 (x,t). 
4J = 1 
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In the following, we will estimate the L p -norm of each term respectively. 
(6.2) 

II^'(''*)IIl p (b p ) 

dx 



< M p 



< M p 



\v\<v 
1 

dx I dz 



\y\ z v(dy) / dz d^ x] v(x + zy,t) 



\v\<v 



dy\y 



= M>[\S 1 (0)\ 



< MP |5i(0)| 



= MP 1^(0) 




2 — 7i — a 



2 — n— a 




y\<v 



v(dy)\y\ 2 \d 2 xZxJ v(x + zy,t)\ 



d X i x M% + zy,t)\ 



y\<v 



dy\y\ 



dx \ d x i x] v(x + zy,t)\ 



y\<v 



dy\yr n - a m*M-,mi p(Bp+v) 



Here the first inequality follows from Fubini's Theorem and Jensen's inequality with respect to the Lebesgue measure 
dz. Assumption (H5) is used in the second inequality. The third inequality follows from Holder inequality with 
l/p+ 1/q = 1. In the second equality, |Si(0)| is the surface area of the unit ball in R". Note that x + zy G B p + n 
when x £ B p , z 6 (0, 1) and \y\ < rj, the fourth inequality follows. 

For I 2 and I3, noting that x + y G B p+1 when x £ B p and \y\ < 1, we have 



(6.3) 
(6.4) 



|/a(-,i 



IL P (B P 



< C ■ max Bp+lX[ ^ T] \ V x v\ ■ 



(I + t? 1 -"), a^l 
(1 - \ogr]), a = l 



and 



ll J 3(-,*)|| jLp(Sp ) < C-max R „ x[aiT] |v| • J |y|>1 v{dy). 



Combining (6.2) - (6.4), (6.1) follows from \\Iv\\ Lp ( B pX(s,T)) = Is \\Iv(;t)\\L p (Bp)dt 
\\d* ixj v\\L p (B p+n x{ s ,T)) < \H w ^(Bp + ,x(^T)) ( sec Definition 4.1). 



and 



□ 



In (6.1), when a £ [0, 1) (finite variation jumps), the factors of r\ in both terms on the right-hand-side converge to 
as r\ — > 0. Therefore, the L p -norm of Iv on the domain B p (xq) x (s, T) essentially only depends on max r x r s _ T ] \v\ 
and max Bp+lX [ Sj j'] |V x w|. This can be also confirmed by working with the reduced integral form I^v in (3.2). 

On the contrary, when a £ [1,2) (infinite variation jumps) , the factor 1 + rj 1 ^ 01 (or 1 — log rf) in (6.1) will blow 
up as 77 — > (a similar phenomenon was also observed in Lemma 1.1 of [6] pp.206 for L p -norm on R"). Therefore, 
it is important to note that the L p -norm of Iv on the domain B p (xq) x (s,T) actually depends on W^-norm of 
v on a larger domain B p+ri (xo) x (s,T). Because of the expansion of the domain, instead of using the boundary 
estimate in Theorem 9.1 in [22] pp. 342, we will use the interior estimation technique in Theorem 10.1 in [22] pp. 
351 to prove Theorem 5.1 in the following. 

Proof of Theorem 5.1. Let us choose a cut-off function ( s (x,t) such that 



(\x,t) 



1 (x,t) eB p x (S,T) 

{x,t)£M n x(0,T)\B p+1 x(^,T) 
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Here the constant S € (0,s) will be determined later. This cut-off function can be chosen such that 



(6.5) 



Co 



c. 



d*'C\<-f> |^C a |<^ and \d t C\<f 



f° r hj< n and some constants C-y, C2 and C3. Please see Figure 1 for the domains used in this proof. 
Defining u(x,t) = £ s (x,t)v(x, t), it satisfies 

(fit ~C D +r) u(x, t) = C 5 ■ Iv(x, t) + C S ■ f(x, t) + h(x, t), (x, t) G B p+ s_ x (0, T), 

u(x, t) = 0, (a;, t) G dB p+s x (0, T), 



u(x, 0) =0, x G B p+ s . 



in which h(x, t) = d t ( 5 ■ v — Y^i j=i a u (^xjC 5 ' v + 2 fi^C 5 • d x jv) — Y17=i ^ ' ^e*C* ' w - Appealing to Theorem 9.1 
in [22] pp.341, we can find a constant C such that 



ll u llw#- 1 (B ii>+ |X(o,T)) - C 



(6.6) 



I I I Jj„ II II II II 1j c 



i.j = l 



in which all _L p -norms on the right-hand-side are on B p+ s x (0,T). 

In the following, we will estimate the terms on the right-hand-side of (6.6) respectively. 

(6.7) 

IK ' Iv \\l p (B p+ sx(0,T)) - W Iv h p (B^ A x(lT)) 



< C 



p+i 
2-a 



\ V \\w^(B p+ sx(IX)) + C [ l 



max If I + max |V x d| 

K"x[0,T] B p+ s +l x[0,T] 



Here the first inequality follows from the choice of the cut-off function £ , the second inequality follows from 

_ 6 

4 " 2 ' 



Lemma 6.1 for a 7^ 1 case by picking 77 = f and s = I. When a = 1, we also have an estimate similar to (6.7). On 



Figure 1. Domains used in this proof 
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the other hand, we have 
(6.8) 

Moreover, we obtain from (6.5) that 



^ ' f\\L p (B p+ sx(0,T)) - H/Hl p (b p+ |x(|,t)) ■ 



(6.9) 



Similarly, thanks to (H7'), we also have 



< max \v\ 



x[0,T} \J B , x(»,T)\B p x(6,T) 

\ p+ 4 
1— p 

< C max \v\ ■ S p . 

R"x[0,T] 



dtdx V 1 
Sp 



V" a d 2 C 5 



L P (B p+A x(0,T)) 



< C max R » x[0jX] \v\ ■ S~ 



L P (B sx(0,T)) 



< C max B xM \ V x v\ ■ S p and 



(6.10) 
(6.11) 

(6.12) UEJUV^C 4 ' v |l Z, p (.B . x(0,T)) — ^ max R n x[0,T] M ' ^ p • 

Plugging (6.7) - (6.12) into (6.6) and noticing the choice of the cut-off function ( s , we obtain 



(6.13) 

liop' 1 (B p x(S,T)) 



p (S p+ |X(0,T)) 



l t 'llw, a - 1 (B p+j x(f 1 r)) + c 



l + <5 



+ \\fh p (B p+1 x(4,T)) • 

Multiplying <5 2 on both hand side of (6.13) and defining 
K{5) = C 



l—p l-2p 

1 - a + 5— +6— 



max It? I + max IV^f 
i"x[0,T] B p+1 + 1 x[o.T] 



max lul + max IV^I 

R»X[0,T] S p+ | + 1 x[0,T] 



52 II/Hl p (b A x(f,T)) . 



we obtain 
(6.14) 



S 2 \\v\ 



w P ^(B p x{S,T)) 



<4C - 



(i) 2 " a -(0 2 |HI^(B +i x ( |, r) ) + ^)- 



Let F(t) = t l' L 'll li , 2 ' 1 ( J B +s x(tT))- ^ ne inequality (6.14) gives us the following recursive inequality 



(6.15) 



F(S)<4C[-) F(S/2) + K(S). 



Since a < 2, we can choose sufficiently small 6 such that 4C (<5/4) 2 Q < |. Therefore, we have from (6.15) that 



(6.16) 



f(a)<-F(*/2) + 
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On the other hand, thanks to the assumption v € W p ' loc (M. n x (0,T)), F(S) is finite for any 6 G (0, So)- Iterating 
the recursive inequality (6.16) gives us 



oo 1 /AX °° 1 

F(6)<Y /Ti K( Ti )<j: Ti K(S) = 2K(S), 



where the second inequality follows from noticing that K(S) is increasing in S. Therefore, it follows from the 
definitions of F(S) and K(S) that 



MIm/Mz-d „,„ < ||«IIt; 



Wp' (BpX(s,T)) 



Wp' L (B p x(S,T)) 



< 2C 



<C S 



1—V 1— 2p 



1 + 5 L - a +5— 



max v 



max iVa-wl 



K»X[0,T] S p+ | + 1 x[0,T] 



\L p (B p+l x(i,T)) 



max \v\ + max I X7 x v\ + II f II , m ^/s^ 
:»x[0,t]' 1 B p+s+l x[o,T] 1 1 IU IIL p( B P+ |X(2^)) 



□ 



Appendix A. Proofs of several lemmas in Sections 2, 3 and 4 



Proof of Lemma 2.1. Throughout this proof, in order to distinguish the Euclidean norm in K" from the absolute 
value in R, we denote the Euclidean norm as || • || and the absolute value as | • |. Actually, the norm || • || is equivalent 
to the sum of the norms | • | among all components, i.e., 



(A-l) 



V 



<^2\y l \ <n\\y\\, for any ye 



Thanks to (A-l), (2.4) and (2.5) can be proved under a slightly weaker assumption (H2) than JLi^ \y\ 2 v idy), 
which is the main assumption of Lemma 3.1 in [27]. We will only prove the second and third estimates in (2.5) in 
what follows. 

Following from (1.1) and (2.2), we have for any r £ 7o,t that 



(A-2) 



\Xf -x\\ < 



f b(X*,s)ds 
Jo 



a{X x s ,s)dW s 



lim J T e 

ej,0 



The difference of our proof from the proof of Lemma 3.1 in [27] is the estimation of the large jump term \\^\ 
We will focus on the estimation of this term in what follows. 
First, it follows from (2.2) and the triangle inequality that 



(A-3) 



E \\jr\\ =E 



J\\y\\>l 



y^(ds,dy) 











/ / yn{ds,dy) 


+ E 




JO J\\y\\>l 





yv{dy) 



all>i 
/ 



Let us estimate the two terms on the right-hand-side of (A-3) separately. On the one hand, J Q J^ y ^ >1 y /j,(ds,dy) 

is a martingale because of (H2). Hence J Jiui^i D Jj>(ds, dy) is a submartingale (see e.g. Problem 3.7 in [21] pp. 
13). It follows from the Optional Sampling Theorem that 



(A-4) 



E 



y fi(ds, dy) 



J\\v\\>i 







< E 


/'./ 




Jo J\ 



yn(ds,dy) 
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Thanks to (A-l), we can estimate the right-hand-side of (A-4) as follows: 



E 



Jo J\\v\\>i 



y^(ds,dy) 



(A-5) 



<^E 

n 
i=l 



t 

J ||; 
/ 



y l n(ds, dy) 



y l n(ds,dy) 



J\\y\\>l 



E/ ds 

t^i J o J\\v\\>i 



y l \ "(dy) 



< 



E f [ f>*| Kds,dy)+ f ds f E|^| u{dy) 
Jo J\\y\\>l i= i Jo J\\y\\>i i= i 

2 f ds ( V|y*| v{dy)<2n( \\y\\ v{ 
Jo J\\v\\>i i=1 «/||w||>i 



'{dy) ■ t. 



Here the first and fourth inequalities follow from (A-l). Moreover, the third inequality follows since the Poisson 
random measure is a non-negative measure on R + x W 1 for each to £ fl. On the other hand, the second term on 
the right-hand-side of (A-3) can be estimated similarly using (A-l). 

Thanks to (A-3) - (A-5), we can find a positive constant C such that E || < Ct for any r € 7o,t- The other 
three terms on the right-hand-side of (A-2) can be estimated in the same way as in Lemma 3.1 of [27]. In particular, 
the stochastic integral and the small jump terms are bounded by C t 1 ^ 2 . Moreover, compared to the estimate (3.3) 
in [27], the boundness of b and a ensures that the constant C in (2.5) is independent of x. 

To prove the third estimate in (2.5), we will still focus on the large jump term. Instead of applying the Doob's 
inequality as in Lemma 3.1 in [27], we will use properties of fi to derive the following estimate: 













E 


sup J"/ 


= E 


sup 






.0<s<t 




0<s<t 





yfi(du,dy) 



< E 



sup 2 

0<s<* .-, 

. — — 2 — 1 



(A-6) 



< E 



sup 

0<s<tJo ^||y||>i" 



E W\ V(du,dy) 







< E 


/'/ 




Jo J\\ 



y l [i{du, dy) 
J\\y\\>l 

n 

E Wlvidv^dy) 



a% a 71 r. 

du E \y l \ v ( d y) ^ n 

Jo J\\y\\>i l=1 J\\y\\>i 



y\\v(dy)-t. 



Here the first and fourth inequalities follow from (A-l), the second and the third inequalities hold since fi is a 
non-negative measure for each w £ fl. The rest proof follows from the same approach used in Lemma 3.1 of 
[27]. □ 



Proof of Lemma 3.1. Wc shall show the Holder continuity in x first. Let us break up the integral into two parts: 

(A- 7) \l? u{x\,t) — I^u(x 2 , t)\ < J K \u(xi + y,t) — u{x\,t) — u(x2 + y,t) + u(x 2 , t)\ v{dy) < I\ + I 2 , in which 
(A-8) h = I\y\< e [\ u ( x i +U,t) -u(xx,t)\ + \u(x 2 + y,t) -u(x 2 ,t)\] v{dy), 

(A-9) I 2 = Jj y | >e [ \u{xi +y,t) - u(x 2 + y,t)\ + \u(xi,t) - u(x 2 ,t)\) v{dy). 

Here the constant e 6 (0, 1] will be determined later. Since x — > u(x,t) is globally Lipschitz (see Lemma 2.2), 

\u(xi + y,t) - u(xi,t)\ < L x \y\, \u(xi + y,t) — u(x 2 + y,t)\ < L x \xi-x 2 \ and \u(xi,t) - u(x 2 ,t)\ < L x \xi-x 2 \, 
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for i = 1,2. Combining these inequalities with (H5), in which < a < 1, we obtain from (A-8) and (A-9) that 
(A-10Xi < / 2L x \y\ v{dy) < 2L X M f \y\^ n - a dy = 2L X M\8 1 (0)\ \\- a dr = 2L *M\Si(0)\ ^ 

J\y\<e V |y|<e Jo 1 _ a 

(A-ll)/ 2 < / 2L x \x x ~x 2 \v{dy) < 2L x \x x - x 2 \ f v{dy) + 2L x M\x\ — x 2 \ \ \y\~ n ~ a dy 
= 2L x \ Xl -x 2 \ v(dy)+2L x M\Si(0)\\x 1 -x 2 \ 

J\v\>l 



>\y\>e J\y\>l J £<\y\<l 

± if < a< 1 
l\ y \>i ' [ -loge if a = 0, 

where |Si(0)| is the surface area of a unit ball in M n . Now picking e = \xi — x 2 \ A 1 and noticing that < a < 1, 
we haw 



(A-12) e 1 -" < \x x -x 2 \ l - a , e- a -l<\x x -x 2 

Moreover, when e = \x\ — x 2 \ < 1, 

(A-13) -loge= / -dz< f ^ d z = ~ (\xi - x 2 \~P - l) < ~\ Xl - x 2 \~P V/3>0. 

Hence choosing e=\x\ — x 2 \ A 1, we have — loge < \\x\ — x 2 \~P for any (3 > 0. Combining (A-7) and (A-10) 
(A-13), we conclude that 



when < a < 1 : \l*u(xi, t) — I^u(x 2 ,t)\ < 



2L x M\Si{ 



2L x d a f v{dy) 

J\v\>l 



xi - X2I 1 Q , 



when a — : \l^u(xx,t) — I-'u(x 2 ,t)\ < 
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a(l-a) 7| y |>i 

B M|5i(0)| + 2L * M f l(0)l + 2L x d^ f v(dy) 
P J\y\>l 



\x\ - x^ 1 ' 



in which (3 £ (0, 1) and d = max^^gn \x — y\. 

Similarly, in order to show the Holder continuity in t, we also break up the integral term into two parts: 

(A-14|/^m(x, t\) — I^u(x, t 2 )\ < J K \u(x + y, ti) — u(x, ti) — u(x + y, t 2 ) + u(x, t 2 )\ v(dy) < I\ + I 2l in which 
(A-15) h = f lyl < e [\u{x + y,t 1 ) -u(x,h)\ + \u(x + y,t 2 ) -u(x,t 2 )\] v{dy), 

(A-16) I 2 = J M>e [ \u(x + y,ti) - u(x + y,t 2 )\ + \u(x,ti) - u(x,t 2 )\ } v{dy). 

The constant e £ (0, 1] will be determined later. We can first bound I\ in (A-15) using (A-10). Then it follows 
from the semi-Holder continuity of t — > u(x,t) (see Lemma 2.2) that 

h < I 2 L t \h - t 2 \* v(dy) = 2 L t \h - t 2 \* [ u{dy) + 2 L t \h - t 2 \?[ v{dy) 

J\y\>e Je<\y\<l J \y\>l 

(A-17) ( -a, , 

if0<a<l 



loge, if a = 0, 



<2L t |ti-t 2 |* / i/(dy)+2L t M|5 1 (0)||ti-t 3 |'- 

" / |y|>i 

in which the second inequality follows from (H5) with < a < 1. 

Now picking e = jti — £ 2 |^ A 1, we have e 1 "" < \ti — t 2 \^~ and e~ Q — 1 < \ti — t2| - ^. A calculation in (A-13) 
gives us that - loge < 2\h -t 2 \- f3 / 2 /l3 for any f3 > 0. Therefore (3.5) and (3.6) follow from combining (A-14), 
(A-10) and (A-17). □ 

Proof of Lemma 3.2. Let us first consider the nonlocal boundary value problem: 

(-d t - C + r) v(x, t) = 0, (x, t)£Bx [ti,t 2 ), 

(A-18) 

v(x,t) = u(x,t), (x,t) £ M" x [t 1: t 2 ] \ B x [ti,t 2 ), 
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in which B is the same domain as in (3.7). The viscosity solution of (A-18) is defined as follows. (See e.g. Definition 
12.1 in [9].) 

Definition A-l. Any v G C°(B x [h,t 2 ]) is a viscosity subsolution of (A-18) if 
(A-19) {-d t -C + r)(f>(x,t) <0, for (x,t) G B x [h,t 2 ), 

(A-20) mm{(-d t - C + r)<p(x,t), v(x,t) ~ u(x,t)} < 0, for (x,t) G dB x [t 1 ,t 2 )UB x t 2 , 

(A-21) v(x,t) < u(x,t), for (x,t) eR n x[ti,t 2 ]\B x[t 1 ,t 2 ], 

for any function <j) G C 2,1 (R" x [ii,^]) H Ci(R n x [£i>*2]) smc/i iftaf <f>{x,t) = v(x,t) and <j>(x,i) > v(x,t) for any 
(x,t) G R™ x [ii , ^2] ■ The viscosity supersolution is also defined analogously. As usual, a viscosity solution is both 
a subsolution and a supersolution. 

Using this definition, it is easy to check that u is a viscosity solution of (A-18). Now the statement of the 
lemma follows from the equivalence of Definition A-l and Definition 3.1. The proof this equivalence is only a slight 
modification of the arguments in the proof of Lemma 2.1 of [30]. □ 

Proof of Lemma 4.1. For the notational simplicity, the constant C denotes a generic constant in different places 
in this proof. 

1. Let us first estimate maxg- 1/</>|. Following (2.10), for (x,t) G Q s , we have 



(A-22)|J0(a:,t)| < 



v{dy)+ / \<t>{x + y,t) - 4>(x,t)\v{dy) 

y| >i 



[x + y,t)- ct>(x, t)-Y,f M(x, t) 

i=l 

max |0| / v(dy) 
J\v\<i~i h»x[o, s ] J M>1 

< H\\ { 3 [ \v\ f> v{dy) + 2 max |0| f v{dy) 

Q> J\v\<l K"x[0,s] J M>1 



I\ V \<1 K"XlU,Sj J\y\>l 

< C ( max |0| + ||0|g 

\R«X[0,s] Qs 

In the second inequality of (A-22), Zi are some vectors in M™ with \zi — x\ < \y\. Therefore, when x G SI, we 
have x + Zi £ ft 1 . The third inequality follows from the Holder continuity of d x i<j> on Q\, i.e., J27=i l^a^O 2 *)*) — 
d x i<f>(x,t)\ < \\<f>\\^ |y| /3_1 . We apply (H5) to obtain the last inequality. Note that (3 > a, hence f M<1 \y\~ n+!3 - a dy 
is integrable. 

The proof of the Holder continuity of x — > I(f>{x,t) and t — > I(j>{x,t) are similar to the proof in Lemmas 3.1. Let 
us check the Holder continuity in x first. For any xi, x 2 G SI and t G [0, s], breaking up the integral term into three 
parts, we obtain 

(A-23) \I(j)(x u t) -I<t>(x 2 ,t)\ <h+h+h, in which 

h(x,t) = / [\4>(xi +y,t)- 4>{xi,t) - y ■ V x (/)(xi,t)\ + \4>{x 2 +y,t) - 4>{x 2 ,t) - y ■ V x (j){x 2l t)\] v{dy), 
J\y\<* 

h{x,t) = / []0(xi +y,t) - 4>{x 2 +y,t)\ + \<f>(x!,t) - <j>{x 2 ,t)\ + \y\ (V^i,*) - V x 0(x 3 ,f)|] v i d v)i 

Je<\y\<l 



h{x,t)= I 

•%l>i 



(xi +y,t)- (f>(x 2 + y,t)\ + \(j){xr,t) - (f>(x 2 ,t)\]v(dy). 
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Here the constant e < 1 will be determined later. Let us estimate each integral term separately. An estimate 
similar to (A-22) shows that 



(A-24) 



h < 2U\M f \yfv(dy) < 2MU\m f \y\- n+ ^ a dy = CU\\§ * 



Thanks to the Lipschitz continuity of x — > 4>{x,t) and the Holder continuity of x — > d Xi <f>(x,t), we can estimate 7 2 
and I3 as 



(A-25)/ 2 < 

< M 



e < I y I < 1 



2 max \V x <j>\\x 1 -x 2 \ + \\<l>M\v\\xi-X2\ l> -* 
R n x[o,s] y s 



v(dy) 



€<\y\<l 



2 max |V^||x 1 -x 2 | + ||0|g|2/||x 1 -x 2 |' 3 - 1 



\y\- n ~ a dy 



e 1 -" - 1 when 1 < a < 2, 



C max iV^l^-x^e-* -l) + CU\\ { J£\x 1 -x 2 \ 13 - 1 

B n x[0,s] Q s — loge when a = 1. 



(A-26YZ3 2 max |Va;<A| Ixi — X2I / v(dy). 

R«x[0, 8 ]' J\y\>i 



Now pick e = \x\ — x 2 \ 1 ^ 2 A 1. Note that 1 < a < 2, we obtain e* 3- ™ < \x± — x 2 \^~ , e~ a — 1 < |xi — x 2 |~^, 
e 1- " — 1 < |a?i — x 2 | i ^ £ ' and — loge < \\x\ — x 2 |~ 5 for any S > (see (A-13)). Since /3 > 1, we will choose 6 = 
in the following. Concluding from these inequalities and (A-23) - (A- 26), we obtain 



08) \, ,i= 



(A-27) \I<f>(xi,t) - I<Kx a ,t)\ < C a [ max |V^| + \\<f>\m ) \ Xl - x 2 

\K"x[0,s] Qi J 

where Co is a sufficiently large constant independent of x±, x 2 and t. 

For the Holder continuity of t — > I<f>(x,t), since <p S H^'^(Ql), it follows from Definition 4.1 that 

n 

^IQc^fotiJ-Qc^fota)! < -hi** 1 , for x e O andti,t 2 G [0,s]. 

Picking e = [xi — Xa] 3 A 1, an estimation similar to Lemma 3.1 gives us 

(A-28) \I<t>(x,h) - I<Kx,1*)\ <Ca(Lt + \\4>\\^) \h-h\^, 

where Co is a sufficiently large constant independent of x, t\ and t 2 . 

Now the first part of the lemma follows from (A-22), (A-27) and (A-28). 

2. Noting that max Ds \<j>\ < \\<f>\\^ and max tl;t2 g [0 , s] M^iit^Mi < s^W^ (see Definition 4.1), the 
second part of the lemma follows from the same argument which we used in the first part of the proof. □ 

Proof of Lemma 5.3. For any Rq > 0, let us consider the following function 

vn 

w ( x > *) = [/ (N) + C i *1 + *). 

where /(ii) = yq-^ and the positive constant C\ will be determined later. It is clear that f(R) is an increasing 
function on (0, +00) and liniR-^+oo f{R) = +00. On the other hand, |9 x »/(|x|)| < ^S^Jp < 1 for any i < n. 
Moreover, one can also check that lim.|. r |_ > ._|_ o0 \d 2 ixJ /(|x|) | = and lim^i^o |^ i2 .j/(|x|)| = 25ij for any i,j < n. 
Therefore both d x if(\x\) and d\ x jf(\x\) are bounded on M. n . Thanks to these properties, we can find an upper 
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bound for |J/(|x|)| as follows: 



|//(W)| = 



i=l 



<dy) 



(A-29) 



< 



/ v{dy) [ dz(l-z) V |yV'| \9'. 
J\v\<i Jo ., , =1 



' 2 . 

x^xi 



f(\x + zy\)\ 



< C 



\y\ 2 v( d v) 



\v\<i 



< +00, 



|2/I>1 



u(dy)\f(\x + y\)-f(\x\)\ 



|y|>i 



for some sufficiently large constant C > 0. Here the last inequality in (A-29) follows from (2.3) and (H2). 
Now, applying the parabolic intcgro-diffcrcntial operator to w, we obtain 



(d t -C D -I + c) w(x,t) >(d t -C 



D 



f(Ro) 



(/(W) + c x t) 



f(Ro) 



ci- E ««^/(W)-E^^/(W)+ c /(i' 

i,i=i 



i=l 



where the first inequality follows from the assumption that (dt — Co — I + c) v(x, t) > 0. We can choose a suffi- 
ciently large constant C\ independent of Rq such that 



(A-30) 



{d t -C D -I + c) w(x, t) > 0, for (x, t) e R" x [0, T]. 



This is because ^ /(M), 9 X » /(M) an d coefficients ay, &i, c are all bounded, c > and |//(|x|)| is bounded. 

On the other hand, w(x, 0) = 7p^y/(M) + v(x,0) > thanks to the assumption v(x, 0) > 0. Moreover, 
when |a;| = i?o, w(x,t) = jtj^a if(Ro) +C\t) + v(x,t) > m + v(x,t) > due to the assumption v(x, t) > —m. 
Furthermore, when \x\ > Rq, we also have w(x, t) > m + v(x, t) > since f(R) is an increasing function. Therefore, 
we claim that w(x,t) > for (x, t) S Br„ x (0, To]. Indeed, if there are some points (x,t) G Br x (0, To] such 
that w(x,t) < 0, w(x,t) must take its negative minimum at some point (xo,io) G Br x (0, To]. Since w(x,t) > 
for x > i?o, we have w(x ,t ) < w(x,t) for all (x, t) G R" x (0,T]. As a result, we obtain d t w(xo,to) < 



°i Y^=ibid x iw(x ,t ) = and 



- 1 «m a 2 , . 



;(xo,to) > (see e.g. Lemma 1 in [14] pp. 34). Moreover, 



Iw(xo,to) > 0, since w achieves its minimum at (xo,to) and V K w(xo,to) = 0. Therefore, we have 

(d t - £d - I + r) w(x ,t ) < 0, 



which contradicts (A-30). 

Now, for any point (x, t) G R™ x (0, T], taking i? 



-00, we have v(x, t) > since lim^^+oo f(Ro) = +00. □ 



Proof of Lemma 5.4. First, thanks to Lemma 5.2, \v e \ is bounded on R™ x [0,T]. In the following, we will show 
it is bounded uniformly in e. It follows from (5.3) (i) that (dt — Cd — I + r)v e = —p e (v e — g e ) > 0. Note that 
v e (x, 0) = g e (x) > (see (4.5)), the first inequality in the statement follows from Lemma 5.3 directly. On the other 
hand, defining u = K + I — v e , u satisfies 



(A-3I) 



(d t -C D -I + r)u = r(K + 1) + p e (v e - g e ), (x, t) G 



(0,T]. 



It follows from (4.5) and (5.3) (ii) that p e (K + 1 — g e ) = with e < eo < 1. Combining with (A-31) and the mean 
value theorem, we obtain 



(A-32) (d t -C D -I + r)u+p e (K + l-g £ )-p e (v e -g e )= d t - C D - I + r + p c (y) u = r (K + 1) > 0, 
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for some jel. Note that both K + 1 — g c and v e — g t are bounded, p e is bounded in any bounded domain. 
Therefore, we have that r + p e (y) is bounded and nonnegative (see (5.3) (iv)). Applying Lemma 5.3 to u and 
picking c = r + p' e (y), we obtain u(x, t) = K +1 — v e (x, t) > on W 1 x [0, T}. □ 
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